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SECTION  I 


INTRODUCTION 

Composite  tube  has  attracted  attention  as  a possible  specimen 
geometry  because  it  offers  several  advantages  over  straight- sided  coupons. 
One  of  the  advantages  is  that  it  can  provide  data  under  combined  ioading 
conditions.  As  a corollary  to  this. the  effect  of  end  constraints  commonly 
observed  in  off-axis  tests  [1]  can  be  eliminated.  In  case  of  angle-ply 
laminates  a tube  is  not  subject  to  the  so-called  free  e^ge  effects  [2]  which 
cannot  be  avoided  in  coupon  apecim.ens. 

Furdiermore.  the  need  for  testing  under  a combined  state  of  stress 
arises  if  one  wants  to  answer  some  of  the  basic  questions  such  as  the 
symmetry  of  the  elastic  tensor  [ 3].  the  equality  of  the  tension  and  com- 
pression moduli  [4].  and  the  interaction  between  the  transvers^:  and  shear 
stress  components  in  the  matrix/interface- controlled  failure  (Sj. 

In  spite  of  dte  aforementioned  advantages  and  needs,  the  available 
data  from  tubular  specimens  is  radter  minimal,  the  main  reason  being 
the  high  cost  of  fabricating  and  testing  tubes  of  high  quality. 

Another  problem  associated  with  tubular  specimens  of  anisotropic 
material  is  that  die  uniformity  of  the  state  of  stress  produced  depends  on 
the  anisotropy  as  well  as  die  geometry.  Tube  dimensions  i squired  to 
achieve  a uniform  state  of  stress  have  been  studied  both  analytically  and 
theoretically  [6-12].  Aldiough  the  exact  geometry  depends  on  the  material 
properties,  it  has  been  found  that  the  gage  section  should  be  at  least  twice  as  long 
as  the  diameter  and  that  the  wall  thickness>to  diameter  ratio  should  be  less 
than  0.03  to  ensure  fairly  uniform  stress  distribution.  Typical  tube  dimen- 
sions reported  in  the  literature  are  listed  in  Table  1.  The  tubular  specimens 
tested  in  the  present  study  have  the  same  dimensions  as  none  of  these. 
However,  they  satisfy  the  minimum  requirements  aforementioned. 
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The  preaertt  report  present*  resulte  from  the  combined  loediog  tests 
of  off-axis,  unidirectional  composite  tubes.  The  dam  are  then  analysed  in 
such  a way  as  to  answer  the  questions  raised  above  regarding  the  material 
properties.  It  is  hoped  that  the  amount  of  data  gathered  is  sufficient  to  infer 
statistically  meaulugful  conclusions. 


SECTION  III 


EXPERIMENTAL  PROCEDURE 


1.  SPECIMENS 

Sp*citn«na  are  unidirectiotaal  S-'ply  graphite /epuxy  (T 300/5208)  tubes 
30.  S cm  lour  and  4.06  cm  in  oa^:eide  diameter.  Each  tube  is  identified  by 
the  iib^'^  bi-ientation  angle  follcmred  by  a specimen  number.  Thus  the 
;<pecimen  -45-2.  for  example,  is  one  of  the  tubes.  The  reference 

coordinate  systems  used  to  define  the  fVoer  orientation  angle  0 is  shown  in 
Figure  1,  where  x i«  parallel  ‘to  the  tube  axis  and  1 is  the  fiber  direction. 

The  tubes  wore  purchased  froiri  the  Whittaker  Corp.  and  had  been  kept  in 
die  room  environment  until  the  test.  This  wafting  period  was  long  enough  to 
allow  equilibrium  moisture  ccntent  in  the  tubes. 

The  fiber  volume  content  was  found  to  be  58  A 55i  from  the  photomicro- 
graphs of  a [0]._  tube.  Examination  of  the  photomicrographs  also  revealed 
oT 

that  the  material  had  unusually  large  void  content  ranging  up  to  2. 1%.  The 
voids  manifest  themselves  in  the  low  transverse  strength,  aa  discussed 
later.  Figure  2 is  a photomicrograph  showing  a void. 

2.  TEST  PROCEDURE 

Six  pairs  of  end  fixtures  were  made  to  grip  the  tubes,  A tube  with  a 
grip  attached  at  each  end  is  shown  in  Figure  3.  The  grip  essentially  consists 
of  two  concentric  cylinders.  Figure  4.  The  gap  between  these  two  cylinders 
is  filled  with  an  adhesive  material  and  then  a tube  specimen  is  slipped  in. 

The  adhesive  is  a mixture  of  Epon  828  (6  parts)  and  Versamid  (4  parts). 

The  complete  setup  was  cured  for  c«ne  hour  at  93*^0. 

All  tests  were  performed  on  an  MTS  closed  loop  testing  system  which 
is  fully  computerized.  The  loading  rate  ranged  approximately  from  0.5  to 
5 MPa/s  and  strains  were  measured  at  the  middle  of  the  gage  section  using 
the  three-element  (0/45/90)  strain  rosettes  (Micro- Measurements  Type 
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EA-06-12SRD-350).  T’ne  strasis-atrain  measuremer.ta  were  Uken  >C 
equal  intervals  up  to  the  maximum  applicid  stresses  and  were  stored  in  the 
computer  memory.  The  data  were  analysed  immediately  after  each  tesl.. 
The  following  procedure  was  adopted  to  calculate  the  elastic  complianct.'^s. 

Suppose  « and  or  are  the  strain  and  stress  component,  respectivelly, 
of  interest.  Then  the  paired  data  i = 1.  2,  •••,  20^  were  fit 

by  a linear  equation  of  the  form 

e^*^  = + d 


The  average  slope  S,  which  is  the  elastic  compliance,  was  then  printed  out 
on  the  printer. 

In  all,  four  different  types  of  tests  were  performed  to  characterise  the 
elastic  properties: 


1.  /ocial  loading  - tension  (a  >0)  and  compression  (ff  <0) 

X ^ 

2.  Torsional  loading  - positive  (a  >0)  and  negative  <0) 

xy  xy 

3.  Positive  combined  loading  (o*  /or  >0)  - (v  >0,v  >0)  and 

3C  3C 

(ff  <0,  (T  <0) 

X xy 

4.  Negative  combined  loading  (v  /o-  <0)  - (o-^>0,  a <0)  and 

xy  X X 

(o-  <0,  <r  >0) 

X xy 
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SECTION  III 


RESULTS  AND  ANALYSIS 


1.  ELASTIC  COMPLIANCES 


Elastic  compliances  measured  from  the  tubes  are  summarized  in 
Table  2.  In  all  cases  the  same  tension- compression  and  pure  torsion  in 
both  directions  were  repeated  three  times,  so  that  the  number  of  measure- 
ments analyzed  is  24  for  each  specimen.  It  should  be  noted  that  one  test 
consists  of  loading  and  unloading,  thus  providing  two  measurements  of  the 
same  compliance.  The  compliances  are  defined  by  the  following  equations: 


e = 

S',  , (T 

+ S',, 

X 

11  X 

16 

e = 

S'  , (T 

+ S' 

y 

21  X 

26 

e = 
xy 

61  X 

*®66 

16  xy  ’ 


(1) 
(2) 
(3) 

Here,  the  subscripts  x and  y are  the  reference  coordinates  for  loading, 
Figure  1. 

For  0-degree  tubes  and  should  vanish  since  the  mater- 

ial is  orthotropic.  The  measurements  are  not  exactly  as  predicted  theore- 
tically; however,  the  data  show  large  scatter,  indicating  that  the  nonzero 
values  are  probably  a result  of  experimental  error  as  well  as  of  the  devia- 
tion from  the  assumed  uniform  state  of  stress. 


Unusually  high  coefficient  of  variation  (C.  V.)  in  S'  and  S'  of  speci- 
ment  15-3  is  due  to  the  large  difference  between  tension  and  compression 
moduli  which  is  in  turn  believed  to  result  from  the  poor  alignment. 

In  order  to  check  the  equality  of  compliance  under  loadings  in  opposite 
directions,  e.g.  tension  and  compression,  S.'i'^)  are  plotted  against  S.'(~)  in 
Figures  5 and  6.  The  superscripts  (+)  and  (-)  denote  the  compliances 
obtained  unt'er  positive  and  negative  loadings,  respectively.  In  these  fig- 
ures the  solid  lines  are  the  linear  least  squares  fit  of  the  data.  That  is. 
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the  data  were  fit  by  an  equation  of  Uie  form 


S!(+)  = aS:(-)  +b  , 


(4) 


and  the  corresponding  parameters  a and  b are  listed  in  Table  3,  together 
with  the  coefficient  of  correlation  r.  In  plotting  the  data  absolute  values 
were  used  when  and  S'/.'^have  the  same  sign;  otherwise,  the  actual 
measurements,  including  sign,  were  plotted,  e.g.  Figure  5(c). 

Figures  5 show  that  the  compliances  measured  in  simple  tension  and 
compression  are  equal  to  each  other  except  for  S'.  I S'  | is  seen  from 
Figure  5(b)  to  be  slightly  higher  in  tension  than  in  compression.  It  is  inter- 
esting to  note  that  | tends  to  be  larger  in  tension  for  the  off-axis 

angles  |o|  > 45°  and  smaller  in  tension  when  |d|  < 45°. 

Figures  6 similarly  show  that  the  compliances  measured  in  torsion 

are  independent  of  whether  the  torque  is  positive  or  negative.  Here  the 

positive  torque  is  defined  to  be  in  the  same  direction  as  is  the  angle  9. 

Thus,  the  positive  torque  results  in  a positive  torsional  stress  (a  >0)  if 

xy 

0 is  positive  and  in  a negative  torsional  stress  ® negative. 

The  best- fit  line  i.n  Figure  6(b)  indicates  that  is  slightly  lower  than 
However,  the  three  points  far  off  to  the  right  are  from  the  specimen 
which  exhibited  appreciable  misalignment.  Although  how  tliis  misalignment 
affects  the  compliance  is  not  exactly  known,  it  is  suspected  that  the  devia- 
tion is  due  to  the  misalignment.  Thus,  if  we  neglect  these  measurements, 
the  results  will  undoubtedly  improve  the  equality  between  and  S'^ " ) . 

CO 

Average  values  of  S'  are  plotted  against  average  values  of  S'  , in 

o 1 16 

Figure  7 to  check  the  symmetry  requirements.  The  data  are  from  Table 
2.  For  the  material  tested,  the  deviation  from  the  symmetry  is  rather 
small. 

Averaging  of  the  compliance  data  can  be  performed  by  using  the  invar- 
iants [13,14].  The  necessary  invariants  are 
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(5) 


'2 

R 


(s;,+  s.2-  2Si^+s.^,/8 


1/2 


= i:'-®'ii"®22'  J 


1/2 


[(Sij  + S.2.2S.2-S./+  4(S.^-Si/]  /8 


(6) 

(7) 

(8) 


In  the  above  equations  S'  , stands  for  the  average  of  S'  , and  S'  in  Table  2. 

16  1661 


The  calculated  invariants  are  shown  in  Figure  8'.  The  mean  values  of 
the  invariants  and  the  corresponding  coefficients  of  variation  are  listed  in 
Table  4.  The  invariants  I and  R-  show  higher  scatter  than  do  the  invar- 

1 M 

iants  I2  and  R^.  In  the  calculation  it  was  assumed  that 


8^2(6)  = (90  - 0) 


because  S^2  measured,  and  that  S^^ 

specimens. 


0 for  0-  and  90-deg. 


The  average  invariants  are  then  used  to  calculate  the  compliances 
through  the  equations 


11 

= 'l"2-*l- 

«2 

(9) 

'22 

- 7,  - 

^2 

(10) 

12 

= I j - I2  + R2 

(11) 

'66 

= 412+4^2  8 

(12) 

where  an  over  bar  denotes  average.  The  resulting  compliances  are  also 
listed  in  Table  4. 

Comparison  between  the  predictions  from  the  average  compliances 
and  the  off-axis  data  is  shown  in  Figures  9.  Aside  from  the  experimental 
scatter,  the  curves  based  on  the  average  compliances  are  in  fairly  good 
agreement  with  the  data. 
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2. 


ELASTIC  BEHAVIOR  UNDER  COMBINED  LOADINGS 


The  Etress  ratios  employed  in  the  combined  loading  tests  are  listed 
In  Table  5.  These  tests  can  serve  as  a check  on  how  sufficiently  the  average 
compliances  describe  the  elastic  behavior.  The  loading  paths  in  the  <r^- 
plane  have  been  described  in  Section  II.  The  resulting  strain- to- stress 
ratios  are  then  analytically  determined  from 


e /(T  = 

X X 

S'  + S'  ,0-  /<T 

11  16  xy  X 

(13) 

A 

It 

S'  + S'  IT  /ff 

IZ  Z6  xy  X 

(14) 

e /c  - 
xy  X 

S',  + S'  cr  /(T 

16  66  xy  X 

(15) 

e /cr  = 

X xy 

S'  or  /(T  + S'  , 

11  X xy  16 

(16) 

e /(T  = 

y xy 

S'  or  /a-  + S' 

IZ  X xy  Z6 

(17) 

e /cr  = 
xy  xy 

S'  , 0-  /cr  + S' , 

16  X xy  66 

(18) 

In  Figures  10  the  calculated  values  are  on  the  abscissa  and  the 
measured  ones  on  the  ordinate.  The  straight  lines  in  the  figures  represent 
a perfect  correlation.  In  all  cases,  the  data  are  scattered  closely  around 
these  lines,  indicating  a good  correlation. 

3.  STRENGTH 

The  stress  components  along  the  material  symmetry  axes  at  failure 
are  listed  in  Table  C.  As  mentioned  before,  these  stresses  are  introduced 
by  applying  the  axial  stress  and  the  torsional  stress  simultaneously 
while  keeping  the  ratio  constant. 

If  any  coupling  between  the  fiber  failure  and  the  matrix/interface 
failure  is  neglected,  then  the  matrix/interface  failure  criterion  can  be 
written  as 
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= Vz  * *'22'2  + *'66'"6  " * O’* 

Thf;  strength  tensor  components  F^,  and  are  then  determined  by 

the  linear  least  squares  method: 

= ["fii)  (20) 

where 


L''66J 


and  jl}  is  the  1 x n column  matrix  whose  elements  are  all  unity.  The  super- 
script (i)  stands  for  the  i- th  data  set.  Note  that  the  total  number  of  measure- 
ments, n,  is  26  in  the  present  case.  The  results  are 

F^  = 3.  376  X lO’^' (MPa)"^ 

F22  = 4.721  X lo"^  (MPa)‘^ 

F^^  = 2.  384  X 10“^  (MPa)"^ 

The  corresponding  failure  surface  in  the  o-^  - (r^  plane  is  shown  in  Figure 

11. 

Table  6 lists  the  value  of  f calculated  for  each  tube  from  the  above 
F's  and  the  failure  stresses.  Since  the  minimum  value  of  f is  less  than 
zero,  i.e. 
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f . 
min 


- 0.6035 


(23) 


the  scatter  of  f i>:  fit  by  a Weibull  distribution  of  the  following  form: 

R = “’’[■("^)  ] 

R is  the  probability  of  the  failure  function  greater  than  f.  The  shape  and 
scale  parameters  determined  are 

Q = 3.055  , f = 1.5265 

and  the  coefficient  of  correlation  is  0.9729. 


The  average  value  of  f obtained  from  the  distribution  (24)  is  only  0.76, 
wiiich  is  much  lower  than  the  u.iity  initially  assumed.  This  is  also  apparent 
in  Figure  12  where  there  are  more  data  points  inside  the  failure  surface 
(f  < 1)  than  outside  (f  > 1).  The  reason  is  because  the  least  squares  fit  places 
more  weight  on  the  higher  stresses  through  second  order  terms  in  the  poly- 
nomial. Thus  lower  stresses  have  less  influence  on  the  strength  tensor 
components . 

The  effect  of  the  longitudinal  stress  component  on  the  matrix/interface 
failure  is  studied  in  Figure  1 3 by  plotting  the  value  of  f at  failure  against  a^. 

The  coefficient  of  correlation  for  the  data  is  only  0.0306,  indicating  very 
little  influence  of  on  the  matrix/interface  failure  within  the  range  of 
applied.  Note  that  the  maximum  a ^ is  less  than  30%  of  the  typical  longitudinal 
strength. 

Now  that  the  failure  function  is  known,  one  can  examine  the  margin  of 

safety  involved  in  the  elastic  property  tests.  The  maximum  value  of  f to  be 

reached  in  each  test  can  be  calculated  by  substituting  the  intended  terminal 

stress  components  in  the  failure  function.  The  results  are  listed  in  Tables 

2 and  5.  In  Table  2,  f and  f are  the  maximum  values  of  f in  the  axial  and 

X xy 

torsional  tests,  respectively,  and  the  superscripts  ( +)  and  (-)  in  Table  5 

stand  for  the  sign  of  c /<r  . 

xy  X 
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Comparison  of  the  values  of  f in  the  elastic  tests  with  those  at  failure 
reveals  that  two  specimens,  -15-2  and  60-2,  did  not  fail  in  the  combined 
elastic  tests  although  in  theory  they  should  have  failed;  i.e.,  the  maximum 
f was  exceeded  in  both  cases.  Thus,  it  may  be  concluded  that  the  value  of  f 
at  failure  depe^-rds  on  the  state  of  stress,  which  is  contrary  to  the  assumed 
failure  criterion.  However,  in  all  the  other  cases,  failure  did  not  occur, 
ae  expected,  when  f was  less  than  that  at  failure.  This  indicates  that  the 
discrepancy  in  the  two  cases  mentioned  above  is  probably  due  to  the  variation 
from  specimen  to  specimen  of  the  strength  tensor  components. 

In  all  the  tests,  failure  initiated  in  the  test  section  in  the  form  of  crack- 
ing along  the  fibers.  Typical  failure  modes  are  shown  in  Figure  14  for  every 
off-axis  angle  tested.  Multiple  fracture  of  60-deg  specimens  is  a result  of 
those  specimens  not  being  able  to  sustain  much  torque  after  the  fracture. 


SECTION  IV 


CONCLUSIONS 


Elastic  compliances  and  the  matrix/interface- controlled  failure  sur- 
face of  a unidirectional  graphite /epoxy  composite  (T300/5208)  have  been 
determined  by  testing  off-axis  tubular  specimens.  Invariants  have  been 
used  to  obtain  the  average  compliances  and  the  linear  least  squares  method 
to  determine  the  failure  surface. 

A good  agreement  is  shown  between  the  prediction  and  the  data  for  the 
elastic  behavior  under  combined  loadings.  The  equality  between  the  tension 
and  compression  compliances  and  the  symmetry  of  the  compliance  tensor 
are  established  to  within  the  experimental  scatter. 

The  failure  surface  in  the  (t.  - v,  plane  is  characterized  by  a second- 
order  polynomial  including  a first-order  term  in  <r^.  This  failure  criterion 
agrees  with  the  experimental  observation  that  the  compressive  loading  per- 
pendicular to  the  fibers  can  increase  the  longitudinal  shear  stress  required 
for  failure  (See  [5]  for  Gr/Ep,  [10]  for  B/Ep,  and  [15]  for  Gl/Ep).  This 
is  not  surprising  if  one  notes  that  the  matrix/interface- controlled  failure 
initiates  from  the  inherent  defects,  such  as  voids  and  partial  debond,  which 
are  aligned  fairly  parallel  to  the  fibers. 

Unfortunately,  a large  scatter  is  seen  in  the  strength  data  and  further- 
more the  transverse  strengths,  both  tensile  and  compressive,  calculated 
from  the  strength  tensor  are  lower  than  what  are  reported  in  the  literature. 
This  difference  is  believed  to  be  due  to  the  unusually  high  void  content  de- 
tected in  the  specimens  used.  However,  the  data  still  show  the  common 
failure  characteristics  expected  of  graphite /epoxy  composites. 


IZ 
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TABLE  4.  AVERAGE  INVARIANTS  AND  AVERAGE  COMPLIANCES 


'2 

’'2 

Ave. , 

(TPa)"^ 

26.43 

34.63 

46.22 

5.510 

C.  V. 

, % 

8.96 

5.29 

5.49 

17.50 

(TPa)"^ 

®.2 

(TPa)'^ 

®22 

(TPa)‘^ 

^66 

(TPa)'^ 

9.  33 

-2.69 

101.77 

160. 56 

TABLE  5.  STRESS  RATIOS  IN  COMBINED  LOADINGS 


Spec.  No. 

<r  / <r 

xy  X 

,(4) 

0 - 2 

± 0.2 

0.05 

0.05 

0 - 3 

± 0.2 

0.05 

0.05 

0 - 4 

± 0.4 

0.01 

0.01 

15  - 1 

dfc  0.4 

0.  17 

0.45 

-15-1 

±0.5 

0.24 

0.38 

-15  - 2 

± 0.4 

0.  19 

-0.07 

-15  - 3 

± 0.2 

0.26 

0.04 

15-3* 

± 0.4 

.081 

0. 19 

-30-1 

± 0.8 

0.72 

0.29 

-30-2 

±0.5 

0.  38 

0.09 

-45  - 2 

±1.67 

0.76 

0.  37 

-45  - 3 

± 0.83 

0.43 

0.10 

-60  - 1 

±1.67 

0.77 

0.24 

60  - 2 

± 0.83 

0.48 

0.02 

-60  - 3 

± 0.83 

0.48 

0.02 

90  - 1 

±1.67 

0.  36 

0.  36 

90-2 

± 0.83 

0.  32 

0.  32 

90  - 3 

± 0.83 

0.  32 

0.  32 

90-5 

± 0.83 

0.  32 

0.  32 

90  - 6 

± 0.83 

0.  32 

0.  32 

* Failure  in  negative  loading. 
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TABLE  6.  STRESS  COMPONENTS  AT  FAILURE 


Spec.  No. 


"1 

MPa 


2 

MPa 


''6 

MPa 


Is  <r  /cr  the  same  as 
xy  X 

one  of  those  in  elastic 
tests? 


0 - 1 

375.4 

0 

44.  38 

0.470 

- 

0 - 2 

0 

0 

48.04 

0.  550 

Yes 

0 - 3 

-375.4 

0 

-39.02 

0.  363 

Yesa 

0 - 4 

0 

0 

-48.60 

0.563 

Yes 

15  - 1 

-259.4 

-18.62 

59.50 

0.687 

Yes 

-15-1 

139.2 

9.99 

37.30 

0.716 

Yes 

-15  - 2 

135.1 

-8.03 

3.12 

-0.238 

No 

-15  - 3 

323.2 

-58.68 

-66.20 

0.689 

No 

15  - 3 

11.49 

4.26 

-9.50 

0.  174 

Yes 

-30-1 

118.4 

12.78 

45.26 

0.997 

No 

-30-2 

193.8 

-85. 15 

-17.79 

0.623 

No 

-30-3 

-11,64 

11.61 

6.69 

0.466 

Yes 

-45  - 1 

-13.79 

25.03 

5.62 

1.148 

- 

-45  - 2 

97.22 

-97.22 

0 

1.180 

Yes 

-45  - 3 

155.4 

-85.81 

34.78 

0.867 

No 

-60  - 1 

55.60 

11.05 

51. 34 

1.059 

No 

60  - 2 

86.56 

-34.81 

-64.91 

0.401 

No 

-60  - 3 

116.1 

-78.02 

78.03 

1.691 

No 

90  - 1 

0 

12.82 

0 

0.  510 

Yes 

90-2 

0 

23.67 

0 

1.064 

Yes 

90  - 2''^ 

0 

20.26 

0 

0.878 

- 

90  - 3 

0 

0 

-74.95 

1. 339 

Yes 

90  - 3'  ^ 

0 

0 

-79.36 

1.501 

- 

90  - 4 

0 

8.04 

0 

0.  302 

Yes 

90  - 5 

0 

8.94 

0 

0.  340 

Yes 

90  - 6 

0 

0 

70. 12 

1.172 

Yes 

a.  Influence  of  o-j  is  neglected. 

b.  Retested  after  failure. 
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ALL  DIMENSIONS  IN  cm  UNLESS  OTHERWISE  MENTIONED 
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Figure  9(e).  S' 
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Figure  11.  Failure  Surface  in  <r,  - 9 Plane 
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